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Abstract 

The leading twist theory of nuclear shadowing predicts the shadowing correction to nuclear par- 
ton distributions at small x by connecting it to the leading twist hard diffraction in electron-nucleon 
scattering. The uncertainties of the predictions are related to the shadowing effects resulting from 
the interaction of the hard probe with > 3 nucleons. We argue that the pattern of hard diffrac- 
tion observed at HERA allows one to reduce these uncertainties. We develop a new approach to the 
treatment of these multiple interactions, which is based on the concept of the color fluctuations and 
accounts for the presence of both point-like and hadron-like configurations in the virtual photon 
wave function. Using the developed framework, we update our predictions for the leading twist 
nuclear shadowing in nuclear parton distributions of heavy nuclei at small x. 
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I. INTRODUCTION 



t he g lobal fits that extract nuclear PDFs from 



In this work, we consider the quark and gluon parton distribution functions (PDFs) in 
nuclei at small values of Bjorken x and their reduction as compared with the incoherent sum 
of the nucleon PDFs because of the phenomenon of nuclear shadowing. Most of experimental 
information on nuclear PDFs comes from inclusive deep inelastic scattering (DIS) with 
nuclear targets which measures the nuclear structure function F2^(x,Q^). For x < 0.05, 
F2a{,x,Q'^) < AF2Ar(x, Q^), which is called nuclear shadowing [F2n{,x,Q'^) is the isoscalar 
nucleon structure function and A is the number of nucleonsl. Because of the factorization 

n 

theorem for DIS (for a review, see Ref. [l|]), which relates F2a{x^Q'^) to nuclear parton 
distributions fj/A{x, Q"^) {j is the parton flavor), nuclear shadowing is also present in nuclear 
PDFs, fj/A{x,Q^) < Afj/N{x,Q^) for x < 0.05, where fj/N{x,Q^) is the PDF of the free 
nucleon. This finds evidence in the results o 
various data on hard scattering with nuclei 

Nuclear PDFs at small x play an important role in the phenomenology of hard scatter- 
ing with nuclei. Their knowledge is required for the evaluation and interpretation of hard 
phenomena in proton-nucleus and nucleus-nucleus collisions at Relativistic Heavy Ion Col- 
lider (RHIC) and the Large Hadron Collider (LHC), in real photon-nucleus interactions in 
ult.apenphe.al coUi.ons at the LHC Ull), and in lepton-nucleu= =catteting at the futute 
Electron- Ion Collider (EIC) |12l. Il3l|. In addition, nuclear PDFs at small x are needed for 
the quantitative estimation of the onset of saturation in ultra high energy interactions with 
nuclei, which can be studied at the LHC and the EIC. 

A comparison of the results of the global fits for nuclear PDFs obtained by various 
groups |2HlO| shows significant discrepancies in the predictions for nuclear PDFs at small 
X (uncertainties of individuals fits at small x are also very large [sl, [ol). The main reason 
for this is that the global fits are predominantly based on fixed-target data that do not 
cover the small-x region (by requiring that > 1 GeV^ is sufficient for the applicability 
of the factorization theorem, one limits x > 5 x 10^'^). In addition, the gluon nuclear PDF 
is determined indirectly from the scaling violations using the very limited data. Therefore, 
the extrapolation of the obtained nuclear PDFs to the low values of Bjorken x that will be 
probed at the LHC and the EIC is essentially uncontrolled. An alternative to the global fits 
is provided by the approaches that attempt to predict nuclear shadowing for nuclear PDFs 
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using the high-energy dynamics of the strong interactions. 

We use the so-called leading twist theory of nuclear shadowing [l^. It combines the tech- 
nique used by Gribov to derive nuclear shadowing for the total hadron-deuteron cross section 

n ^ n 

at high energies [lb] and the QCD factorization theorems for inclusive [1] and diffractive [16[ 
DIS. The numerical predictions employ the results of the leading twist QCD analyses of hard 



diffraction in lepton-proton DIS at HERA 17H19|. Although in the leading twist approach 



the hard probe interacts with one parton of the nucleus, in the target rest frame, nuclear 
shadowing appears as the effect of multiple interactions of the projectile (virtual photon) 
with several (all) nucleons of the target. The interaction with N = 2 nucleons is related in 
a model-independent way to the diffractive PDFs of the nucleon. The account of the inter- 
action with > 3 nucleons is model-dependent and sensitive to the underlying dynamics of 



the hard diffraction. The recent HERA data 



17H19| revealed that the energy dependence 



of the hard diffraction in DIS (dependence on the light-cone fraction xjp) is close to that of 
the soft processes. This indicates that the hard diffraction in DIS is dominated by large-size 
hadron-like configurations in the photon wave function. This observation allows us to im- 
prove the treatment of the contribution to nuclear shadowing coming from the interactions 
with A^ > 3 nucleons as compared to the simplified quasi-eikonal approximation used in 
our earlier papers 20|] by taking into account the presence of both point-like and hadron- 
like configurations in the virtual photon. The goal of the present Letter is to present a 
new, improved treatment of such multiple interactions and to update predictions for nuclear 
shadowing in nuclear PDFs. 

We emphasize that the presence of the small-size (point-like) configurations and, in gen- 
eral, configurations of different transverse sizes that interact with different cross sections (we 
call such configurations color or cross section fluctuations) is much more important for the 
virtual photon than for hadronic projectiles. 

Our Letter is organized as follows. In Sec. [Tll we briefly review the leading twist theory of 
nuclear shadowing and present a new formalism for the treatment of the multiple interactions 
of the virtual photon with the nucleons of the nuclear target, which is based on the concept 
of color (cross section) fluctuations. In Sec. Illlt we show that the complete treatment of color 
fluctuations can be well approximated by the so-called color fluctuation approximation. In 
Sec. llVt we present updated predictions for the effect of nuclear shadowing in nuclear parton 
distributions. We focus on the predictions for heavy nuclei at small x since modifications 
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of the predictions for light nuclei (where the N = 2 term dominates) are rather small. Our 
results are summarized in Sec. |Vl 



II. LEADING TWIST THEORY OF NUCLEAR SHADOWING AND CROSS SEC- 
TION FLUCTUATIONS FOR MULTIPLE INTERACTIONS 

The phenomenon of nuclear shadowing is fairly well-understood: in the target rest frame, 
nuclear shadowing arises as the result of multiple interactions of the projectile (virtual 
photon) with several nucleons of the nuclear target. The number of the interactions increases 
with decreasing Bjorken x, which is a result of the space-picture of the strong interactions, 
see e.g., Ref. 2l|. At sufficiently high energies (small Bjorken x), the virtual photon can 
interact with all the nucleons of the target that are located in the photon's path. 

The graphs that contribute to the nuclear structure function F2a{x, Q"^) are presented in 
Fig. [H where we used the optical theorem to relate the imaginary part of the 7* A forward 
scattering amplitude to the nuclear structure function. In the figure, graphs a, b, and 
c correspond to the interaction with one, two, and three nucleons of the nuclear target, 
respectively. Note that the graphs for the interaction with four or more nucleons are not 
shown, but assumed. Graphs b, c and higher scattering terms are responsible for nuclear 
shadowing in F2a{x,Q^). 

The contribution of graph a, which is conveniently denoted F2'^{x,Q'^), is 

Fi'J{x,Q^) = AFM^,Q'), (1) 

where F2n{x,Q'^) is the isospin-averaged structure function of the nucleon. In Eq. ([T]), we 
neglected the deviation from the many-nucleon approximation for the description of nuclei 
and the Fermi motion effect, which are numerically unimportant at small x. 

The calculation of the contribution of graph 6, F2^^(a;, Q^), is fairly straightforward, but 



lengthy 22|. The detailed derivation, including the effect of the real part of the diffractive 
amplitude, is given in Ref. 20|. Here we present the final result for F2^^(a;, Q^), 

H%,Q') = -SnA{A - l)^e^-^—^Bam ^ dxjpF^^^\x,Q\xjp) 

/POO poo 
d% / dzr / dz2PA{b, zi)pA{b. ^^y{^.-^.)-ipm^ ^ (2) 
J —00 J z\ 
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a) b) c) 



FIG. 1: Graphical representation of the nuclear structure function F2a{x,Q^). Graphs a, b and 
c correspond to the interaction with one, two and three nucleons, respectively. The latter two 
graphs and the interaction with four and more nucleons (not shown) lead to nuclear shadowing. 
The dashed vertical lines represent taking of the imaginary part. 

where F^^^"^ is the nucleoli diffractive structure function measured in hard 7*p inclusive 
diffraction; is the nuclear density; rj ^ 0.17 is the ratio of the real to imaginary parts 
of the 7*p diffractive amplitude; -Bjiff = 6 GeV~^ is the slope of the t dependence of the 
diffractive 7*p cross section; xjp is the light-cone fraction of the nucleon momentum carried 
by the Pomeron (see the discussion below); is the nucleon mass. 

The nuclear density pa depends on the transverse coordinate (impact parameter), b, and 
the longitudinal coordinates, zi and Z2, of the interacting nucleons. The ordering Z2 > z\ 
follows from the space-time evolution of the scattering process. The ^A-'^^-^-^)^p~^n factor 
accounts for the excitation of the intermediate diffractive state denoted by X in Fig. [H At 
high energies, the 7*A^ interaction that leads to nuclear shadowing is diffractive in character. 
This is represented by the zigzag lines in Fig. [H It is convenient to think of the zigzag lines 
as depicting effective Pomeron exchanges. In this case, xjp represents the light-cone fraction 
of the nucleon momentum carried by the Pomeron, xjp ~ [M\ + Q^) / + Q^); where Mx 
is the invariant mass of the diffractive state X, and W is the invariant 7*p energy. The lower 
limit of integration over a; ip in Eq. ([2]) corresponds to Mx = 0; the upper limit is determined 
by the typical cut on Mx, Mx < O.IH^^, which arises because of the nuclear form factor. 
Note that Eq. ([2]) is valid independently of the validity of the leading twist approximation 
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a) b) c) 



FIG. 2: Graphical representation of the multiple scattering series for the quark distribution in 
a nucleus. Graphs a, b and c correspond to the interaction with one, two and three nucleons, 
respectively. The latter two graphs and the interaction with four and more nucleons (not shown) 
lead to nuclear shadowing. 



for the hard diffraction. 



One of the key features of the leading twist theory of nuclear shadowing jl4j . |20| is the 
possibility to predict nuclear shadowing at the level of parton distributions. Using the QCD 
factorization theorems for inclusive DIS and hard diffraction in DIS, one can replace the 
observable structure functions by the corresponding parton distribution. This is shown in 
Fig. [21 which represents the multiple scattering series for the quark distribution in nuclei. 
A similar graphical representation can also be given for the gluon distribution using a hard 
probe directly coupled to gluons. 

The contribution of graph a in Fig. [2l which we denote Q^), is readily obtained 

from Eq. ([1]), 

xf^';\{x,Q') = Axf,/^{x,Q'). (3) 
The contribution of graph b is obtained from Eq. ([2]), 

xf%{x,Q^) = -8n AiA - l)^e 'jf i^diff dxjpf3ff^''\/3, Q\ x^) 

/CXD /"OO 
dzr / dz2 PA{b, z,)pA{b, Z2)e'^-''-'^^" , (4) 
-oo J zi 

where f^^^^ is the diffractive parton distribution of flavor j in a proton. According to the 
factorization theorem 16|, ff^^^ is a leading-twist distribution, whose Q"^ evolution is given 



by the DGLAP equations. This is supported by the analyses of diffraction at HERA 
19]. Therefore, the contribution of a;/^y]^(x, Q"^) to nuclear shadowing is also a leading twist 



function, which gives the name to the present approach — the leading twist theory of nuclear 
shadowing. 

The derivation of the expressions for and is general and model-independent: 

the only simplifying approximations are the neglect of nucleon correlations in the nuclear 
wave function and of the t dependence of the elementary diffractive 7*A^ — > XN amplitude. 

Graph b in Figs. [1] and [2] approximates well nuclear shadowing in the low nuclear den- 
sity limit, when the interaction with only two nucleons is important. As one decreases x, 
graph c and higher rescattering terms also become progressively important. To evaluate 
their contribution, one needs to model the interaction of the intermediate state X with the 
nucleons of the target. Our approach is based on the high-energy formalism of cross section 



fluctuations 



23 



26| . which provides a good description of the total hadron-nucleus cross 



sections and, which is far less trivial, of the coherent inelastic diffraction in hadron-nucleus 



scattering; for a review and references, see Ref. 26[. In this formalism, the wave function 
of a fast projectile (virtual photon) is expanded in terms of eigenstates of the scattering 
operator, \a). Each eigenstate interacts with target nucleons with a certain cross section a. 
The probability for the incoming virtual photon to fluctuate in a given eigenstate is given 
by the distribution Pj{<y)- We explicitly show the dependence of Pj{<y) on the parton flavor 
j as a reminder that DIS probes a particular parton distribution of the target. 

The entire series of multiple interactions shown in Fig. [2] can be summed by the stan- 
dard Glauber formalism generalized to include cross section fluctuations, see, e.g., Ref. 26 |. 
Assuming that 1 such that the multiple interactions can be exponentiated, we obtain 

Xfj/N{X,Q^) [ ,2./(. _ -4(l-^„).T^(6)^ 



^/.7A^(^>Q')osr> /■.2,Er=2(-|(l-^^m(&))'^(a'=), 



{a)j J k\ 

where T^(&) = J^^dzp^ib, z); (. . . )j denotes the integration over a with the weight Pj{cr). 
The interaction with k nucleons probes the kth moment of the distribution Pjicr), {cr^)j = 
j-daP,{a)aK 

Equation ([5]) is valid at small x (high energies), when the effect of the finite coherence 
length (the coherence length is proportional to the lifetime of the fluctuations \a)) is unim- 
portant. In this case, the e'^i^i-^2)mNXp fg^Q^Qj- Eq. (jlj) can be set to unity. 



The contribution to nuclear shadowing from the interaction with N = 2 nucleons is 
related in a model-independent way to the diffractive PDFs of the nucleon. This means that 
((T^)j is proportional to the diffractive parton distribution Q, Q, 

Therefore, Eq. ([5]) can be written as 

/ ('e-4(i-*'')-^^W - 1 + 1(1 - i'q)aTA{h)) ) 

-xf,/Nix, Q')aiix, Q')2 / d%— — ^ . (7) 

J V"^ )j 

In order to cast Eq. ([7]) into the more standard form 201], we reintroduce the dependence 
on the longitudinal coordinates zi and Z2, use the definition of aKxjQ"^) from Eq. (E]), and 
identically rewrite Eq. ([7]) in the following form, 

xfj/Aix, Q^) = Axfj/Nix, Q"^) 

\2 /-O.l 



- STTA^gfJe^^-^Ediff /' dxjpf3f^^'\f3,Q^xjp) 



^2 -^{l~ivWi:^dz'pA{b,z') 



X d b dzi dz2pA{b,zi)pA{b,Z2) . (8) 

J J -OO J Zl \^ )j 

Finally, we restore the effect of the finite coherence length by reintroducing the e^i^i-^i)^p'^N 
factor, replace by A{A — 1) to have the correct number of the nucleon pairs, and obtain 
our general expression for the nuclear parton distribution modified by nuclear shadowing, 

xfj/Aix^Q"^) = Axfj/N{x,Q'^) 

\2 /-O.l 



87rA(A - l)3fte^^-^Sdiff /' dxjpl3ff'^^\l3,Q^,xjp) 



oo 



/ 2-4{i-irihf:['dz'pA(b,z')\ 

dz2PAib, z^)pAib, z2)e'i-^-'^>p^^ 1 __ ^ . (9) 

The evaluation of nuclear shadowing in the leading twist theory of nuclear shadowing, and 
Eq. ^ in particular, does not take into account the possible ultra high-energy branching of 
the diffractive exchange which would couple to different nucleons of the target (the so-called 

n 

triple Pomeron fan or enhanced Reggeon diagrams) [27[ . Using the model for the interaction 
with > 3 nucleons which takes into account such diagrams, nuclear shadowing in nuclear 



PDFs was predicted in Refs. 
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29|. 



The general form of the distribution Pj{cr) that enters Eq. (Q is not known. However, 
one can still infer the properties of Pjicr) that determine the strength of nuclear shadowing. 
For the virtual photon, Pj{cr) is very broad and includes the states \a) that correspond to 



both small and large cross sections cr 
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3l|. The fluctuations with small cross sections 



constitute the perturbative contribution to the photon-nucleon cross section; the fluctuations 
with large cross sections correspond to the hadronic component of the virtual photon. In 
practice, (a^) is dominated by the hadronic-size configurations. This expectation is based 
on the QCD aligned jet model 32|, and agrees well with the final analyses of the HERA data 
on hard diffraction which find that ajp(t = 0) = 1.111 ± 0.007 17|, ll8|], which is practically 
the same as in soft processes, q;^^*(0) = 1.0808 [33]. Hence, the diffractive state X in Figs. [T] 
and |2] is dominated by the large-cr hadron-like fluctuations. 

The key feature of Eq. iQ is that it separates the contributions of the small and large 
cross sections [this was the main purpose of rewriting Eq. ([5]) in the form of Eq. ([7]) which 
led to Eq. (j9])]. While the fluctuations with large cross sections contribute to all moments 
{cr''), the fluctuations with small cross sections contribute significantly only to (cr) and (cr^), 
i.e., to the Axfj/]\f{x, Q'^) term and the double scattering term proportional to f^^^\ There- 
fore, since the (. . ■)j/{o''^)j term in Eq. (|9]) probes the higher moments of Pjic), {a'') / (a^) 
with A; > 3, it can be evaluated with the distribution -Pj(cr), which neglects the small-cr 
perturbative contribution and uses only the information on cross section fluctuations from 
soft hadron-hadron scattering. In particular, we assume that the relevant Pj{(y) is equal to 
the distribution over cross section fluctuations for the pion. 



III. THE COLOR FLUCTUATION APPROXIMATION 

The (. . . )j/ (cr^)j term in Eq. ^ can be identically expanded in terms of {cT'')j/ 
with k > 3. We have just explained that the required distribution Pj{<j) is dominated by 
soft hadron-like fluctuations. For such fluctuations, the dispersion of Pjicr) does not lead 
to significant modifications of higher moments of Pjicr), and it is a good approximation to 
use ~ {{(^^) j / {^'^) jY^'^ foi' ^11 k > 3, which we shall call the color fluctuation 

approximation. Therefore, the (. . . )j/ {o'^)j term in Eq. ([9]) is expressed in terms of a single 
cross section, al{x,Q^), 

aiix,Q') ^ {a'),/ {a'), = {{a'),/ {a'),)'^^'-'^ . (10) 
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Applying the color fluctuation approximation to Eq. (Q, we obtain our final expression 
for the nuclear parton distributions modified by nuclear shadowing, 



xfj/A{x,Q ) = Axfj/N{x,Q 

fj/N{ 

d% I ' dzi r' dz2PA{h,z^)pA{h,Z2y^''-'^^'''''"'''e~^'^^-'^^^^^^ (11) 



- xf,,M{x, Q^)MA{A - 1) gfJe^^^^i?diff f '^ dxjpPff^'\p, Q\ xjp) 



In the treatment of multiple rescatterings in the leading twist theory of nuclear shadowing 



in Ref. 20|, one used the so-called quasi-eikonal approximation, which prescribes the use 
of aljyX.Q'^) = a2{x,Q^) in Eq. f|TT]) . While the quasi-eikonal and color fluctuation approx- 
imations give identical results for the interaction with two nucleons of the nuclear target, 
the color fluctuation approximation provides a more accurate treatment of the interaction 
with three and more nucleons. In particular, the interaction with three nucleons is treated 
exactly in the formalism of cross section (color) fluctuations. As follows from the definition 
and modeling (see below) of the effective rescattering cross section crl{x, Q"^), the color fluc- 
tuation approximation represents the scenario corresponding to the lower limit on nuclear 
shadowing within the framework of the leading twist nuclear shadowing. 

To model the distribution Pjic), we assume that Pj{cr) = Pt^{(t), where P-rrio') is the 
distribution over cross sections for the pion. This assumption attempts to capture the ob- 
servation which we explained above that for the higher rescattering contributions to nuclear 
shadowing, only large-size fluctuations of the virtual photon wave function are important. 
This is a somewhat extreme assumption which results in the smallest nuclear shadowing 
(see Fig. H] below). 



The distribution Pn{cr) is conveniently parameterized in the following form 34l |: 



Pj{a) = P^{a) = Ne . (12) 

The parameters N, ctq and Q are constrained by following requirements: 

daP.„{a) = 1 , 
daP^{a)a = a^^^iW'), 

/"OO 

/ daPM^' = {al^{W')Y{l + u^{W')) , (13) 
Jo 





oo 





fOO 
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where is the total pion-nucleon cross section; is the parameter characterizing the 



dispersion of the distribution Pt,{o). Both al^^ and Wo- depend on W = Q /x — Q + m 



In our numerical analysis, we used the Donnachie-Landshoff parameterization for cr^ 



N 
tot 



0", 



TT + N 

tot 



■K-N 
tot 



13.63 (ly 



2n 0.0808 



+ 31.79 (ly^ 



-0.4525 



mb - 



(14) 



Note that in our calculations, we effectively use only the first term in Eq. ([H]), see also 
Fig. Ei The parameter decreases with increasing energy 35|, which means that cross 
section fluctuations decrease with increasing energy. For the pion projectile, u^- ~ 0.4 at the 
pion energy of E.,, ^ 300 GeV 



25 



3|, which corresponds to ^ 600 GeV^. At the CDF 
energy of = (546)2 ^ 3 x 10^ GeV^o;^ ~ 0.16 x (3/2) = 0.24, where the factor 0.16 is 
Ucr for the proton at the CDF energy 25| and the factor 3/2 reflects the constituent quark 



counting 3J]. Assuming a simple linear interpolation between the two energies, we arrive 
at the following model for Ua'- 



0.4-0.16 



(15) 

T4l) and ( IT5l) fully determine 



where = 600 GeV^ and = 3x10^ GeV^. Equations (Il3]), 
Pvr((j) and aKxjQ"^) and their energy (Bjorken x) dependence. In order to keep track of 
our modeling Pj{cr) = Pn{<^), it is also convenient to introduce the notation a^^°^{W'^) = 
aiix,Q'). 

Figure [3] presents a^^°'^(W'^) = al{x,Q^) and cr2{x,Q^) as functions of Bjorken x at fixed 
Qq = 4 GeV^ {W^ = Ql/x — Ql + m\). The left panel corresponds to the -u-quark; the 
right panel corresponds to gluons. (Note that (ji^{x, Q^) is flavor-independent in our model.) 
The fact that af°^{W'^) > al{x,Ql), which is equivalent to {o'^) j / {c'^) j > (o"^)j/(o")j, is a 
general property of the distribution Pj{(j). 



IV. PREDICTIONS FOR NUCLEAR PDFS 



Equation ffTTj) allows one to calculate nuclear PDFs modified by nuclear shadowing. The 
key inputs for this calculation are the diffractive parton distributions f^^^^ and the slope of 
the diffractive 7*p cross section -Bdifi- The current experimental uncertainties of f^^^^ and 
i?diff lead to an uncertainty in the predictions of nuclear shadowing which is much larger than 
the uncertainty associated with the use of the color fluctuation approximation instead of the 
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FIG. 3: The cross sections crf'°"(W2) = cj^(x,Q^) [Eq. (HU])] and o-^(x,Qg) [Eq. as functions of 
Bjorken x at Qq = A GeV^. The left panel corresponds to the u-quark; the right panel corresponds 
to gluons. 



complete treatment of color fluctuations. It is also important to note that cr^ix, Q^) depends 
weakly on energy (Bjorken x). Hence, measuring nuclear shadowing with one nucleus at e.g., 
X = 10^^, will further improve our predictions for all x and A. 

Figure H] presents the ratio of the nuclear to nucleon parton distributions, 
fj/A^x, Q'^)/[Afj/iy{x, Q^)], as a function of Bjorken x at the input scale Ql = 4 GeV^. The 
solid curves correspond to the color fluctuation approximation, Eq. (ITTl) . and o"3(a;,Q^) = 
(J3^°"(H^^), Eqs. (fT3|) . ( IT4|) and (fT5|) . The dotted curves are obtained using the quasi-eikonal 
approximation with (13(0;, Q^) = crli^^Q'^) in Eq. f lTT]) . The two left panels correspond to 
M-quarks; the two right panels correspond to gluons. Note that we added the effect of an- 
tishadowing for the gluon distribution for 0.03 < x < 0.2 using the method described in 
Ref. 20[. The panels in the top row are for ^"^Ca; the bottom panels are for ^"^^Pb. 

As one can see from Fig. HI the color fluctuation approximation corresponds to the 
smaller nuclear shadowing at the input scale Qq = 4 GeV^ than the quasi-eikonal ap- 
proximation. For instance, at x = 10"'^ and for "^^Ca, UA{x,Ql)\c{/uA{x,Ql)\q^e = 
1.10 and gA{x,Ql)\c{/gA{x,Ql)\qc = l-H (the subscripts indicate the color fluctua- 
tion and quasi-eikonal approximations, respectively). At a; = 10~^ and for ^"^Pb, 
UAix, QDici/uAix, (5g)|qe = 1.31 and gA{x, Ql)\c{/gAix, Ql)\cie = 1-29. 

We described the modifications of the nuclear PDFs at small x at the input scale Qq = 4 
GeV^. Predictions for nuclear PDFs at higher scales > Ql are obtained using the 
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FIG. 4: The ratio of the nuclear to nucleon parton distributions, fj /a{'^^Q'^) /[■^fj/N{x-,Q'^)], as a 
function of Bjorken x at the input scale Qq = 4 GeV^. The solid (dotted) curves correspond to 
the color fluctuation (quasi-eikonal) approximation. The left (right) panels correspond to u-quark 
(gluon) distributions. 



standard DGLAP evolution equations with the input given by Eq. (fTTl) . As one increases 
at a given x, the difference between the predictions of the color fluctuation and quasi- 
eikonal approximations reduces: In the DGLAP evolution, the PDFs at small x and large 
are obtained from the region of larger x and smaller Q^, where nuclear shadowing is 
smaller. 

One should also keep in mind that while we present our predictions for nuclear PDFs all 
the way down to x = 10"^, our results for x ^ 10~^ should be considered only as guiding ones 
since, for this kinematic region, various effects beyond the leading twist DGLAP equation 
should start becoming important. 
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V. SUMMARY 



The leading twist theory of nuclear shadowing allows one to predict nuclear parton dis- 
tributions modified by nuclear shadowing. Nuclear shadowing arises as the effect of the 
interaction of a hard probe with a parton simultaneously belonging to several nucleons of 
the nuclear target. While the leading twist theory of nuclear shadowing gives unambiguous 
predictions for nuclear shadowing in the case when the interaction with only two nucleons is 
important, the interaction with > 3 nucleons requires a more detailed knowledge of the 
dynamics of the hard diffractive processes. We propose a new approach to the treatment 
of such multiple interactions using the formalism of cross section (color) fluctuations which 
allows us to take into account the presence of both point-like and average hadronic-size 
conflgurations in the virtual photon wave function. For practical applications, we propose 
a new approximation — the color fluctuation approximation — which approximates well the 
complete treatment of color fluctuations. This approximation improves the treatment of the 
multiple interactions in the quasi-eikonal approximation, which works well in soft processes, 
where color fluctuations are smaller and have a different structure than in DIS. Using the 
developed framework, we present updated predictions for the effect of nuclear shadowing in 
nuclear parton distributions of heavy nuclei at small x. 
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